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Abstract 

We use analytic and numerical methods to determine the density of states 
of a one-dimensional electron gas coupled to a spatially random quasi-static 
back-scattering potential of long correlation length. Our results provide in- 
sight into the 'pseudogap' phenomenon occurring in underdoped high-T c su- 
perconductors, quasi-one-dimensional organic conductors and liquid metals. 
They demonstrate the important role played by amplitude fluctuations of 
the backscattering potential and by fluctuations in gradients of the poten- 
tial, and confirm the importance of the self-consistency which is a key feature 
of the 'FLEX'-type approximations for the electron Green's function. Our 
results allow an assessment of the merits of different approximations: a pre- 
vious approximate treatment presented by Sadovskii and, we show, justified 
by a WKB approximation gives a reasonably good representation, except for 
a "central peak" anomaly, of our numerically computed densities of states, 
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whereas a previous approximation introduced by Lee, Rice and Anderson is 
not as accurate. 



I. INTRODUCTION 

"Pseudogaps" , i.e. a suppression of the low-energy electronic density of states due to 
an interaction effect, are of current interest in the context of high-T c superconductivity^ 
and of low-dimensional organic materials!! and are relevant to the theory of liquid metalsi, 
and to the issue of spin waves above T c in ferromagnetsi. One general mechanism for 
producing pseudogaps involves the presence of long but not infinite range order, for example 
of the superconducting or density wave type. Consider, for definiteness, a one dimensional 
material exhibiting long ranged charge density wave order at T = 0. At T > thermal 
fluctuations prevent long ranged order but at low T the order parameter fluctuates very 
slowly in space and time. It therefore has significant amplitude to back-scatter electrons. 
The back-scattering will tend to open a gap but will not do so completely because the 
order is not perfect. Very similar issues come up in the physics of liquid metalsi, where 
the ions are highly correlated and move slowly in comparison to the electrons, in the theory 
of superconductivity, where slow Cooper pairing fluctuations mix particle and hole states 
(instead of left and right movers) and again tend to open a gap, and in ferromagnets at 
temperatures slightly above T c , slow magnetization fluctuations may suppress electron-hole 
spin excitations and allow weakly damped spin waves to existS. 

A crucial, but still unresolved, question concerns the proper method of calculation of 
physical quantities in the presence of long-range-correlated scattering. A simple and phys- 
ically appealing approximation was proposed in 1973 by Lee, Rice and Andersonll, who 
argued that one need only consider the leading perturbative one-boson-exchange diagram 
for the self-energy. Within this approximation the low energy density of states was constant, 
but was suppressed from the non-interacting value by a factor proportional to the inverse of 
the correlation length of the scatterers. Subsequent workers used different approximations 
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and obtained different results. In particular, Sadovskii obtained a continued-fraction expres- 
sion which led, for the problem studied ini to a low energy density of states which varied 
as the square root of the correlation lengtlJt However, Sadovskii's results have recently 
been called in to question by TchernyshyovB who exhibited a class of diagrams neglected 
by Sadovskii. A yet different set of approximations, the 'fluctuation exchange' or 'FLEX' 
approximation! has been employed by many authors to study pseudogap effects in models 
of superconductivity^, and these treatments have in turn been questioned00. 

In order to clarify this situation, we present results of a thorough numerical and ana- 
lytical study of a simple model exhibiting pseudogap effects, namely electrons moving in 
one spatial dimension and coupled to a backscattering potential which is constant in time 
but slowly varying in space, with correlation length £. Our main results are that (1) the 
approximate treatment of Sadovskii is a good approximation to the numerically calculated 
density of states in the "high energy" regime E > l/^ x , with x an exponent which de- 
pends on whether the potential is commensurate or incommensurate. The terms omitted 
in Sadovskii's derivation noted ini evidently produce negligible corrections. (2) Neither the 
'FLEX' approximation nor non self-consistent treatments such as that of Lee, Rice and 
Andersoni reproduce at all well the dependence of density of states on correlation length, 
although the self-consistency which is an essential feature of the FLEX approximation leads 
to a better representation of the density of states than does the approximation of Lee, Rice 
and Anderson. (3) Fluctuations in the amplitude of the scattering potential play a crucial 
role in the form of the low energy density of states. Most of the models studied make as- 
sumptions about the amplitude fluctuations which are not physically reasonable. (4) Some 
features of the results are controlled by fluctuations in the gradient of the backscattering 
potential; the widely used Lorentzian form has divergent gradient fluctuations which change 
the form of the results. 

The balance of this paper is organized as follows. In section II we present the models 
and a simple 'WKB' treatment which reveals the essential physics. In section III we present 
the results of our numerical study. In section IV we compare the numerical results to those 



obtained by approximate methods. Section V contains a summary, conclusions, and a list 
of open problems. 



II. MODEL 

We consider spinless electrons moving in one spatial dimension and coupled to a static 
potential with spatial correlations to be specified below. The fundamental Hamiltonian is 

h = - { 2 (4 d j+i + h - c ) + + v A d i } (!) 
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We consider two sub-cases: 

A) Half filling and commensurate potential. In this case \x = and we assume V(j) is a 
real random potential chosen from a distribution which implies 

< V 3 >= (2) 

and 

< V s V j+k >= (3) 

Here f(x) — > for x — > oo and we are interested in the large £ limit. 

B) Incommensurate band filling and potential. In this case we linearize Eq. p] about the 
Fermi points, separate the electrons into right-moving (R) and left-moving (L) branches in 
the usual way, and adopt a matrix notation, defining the 2x2 matrix G via 



G R , R (z, t; z', = e <M*-*0 ^ Cfl ( Z) t)4(zV)) (4) 

G^(z, t; z', = c -*^(— ( T c L (^, t)c£ (^)) (5) 

t; = e**^ (Tc^z, i)4(zY)) (6) 

Glr{ z i f; z', H) = e-*'^) ( Tc t t)c R (zY)) (7) 

The Green function obeys a Schroedinger equation which after Fourier transformation in 
time may be written: 



[u + ia 3 v F d z + V(z)\ G = 5(z — z) 



(8) 



A general potential V{z) has a forward scattering part which couples right-movers to 
right-movers and left-movers to left-movers, and a back-scattering part which couples left 
movers to right movers. The forward scattering part is not important for the density of 
states^; thus we consider only the back-scattering part and write 



Viz) = vAz)ai + v 2 (z)a 2 



(9) 



Because of the phase factors arising from the incommensurate Fermi surface V is a complex 
potential; further we have absorbed phase factors into the definition of G, thus we write the 
V correlator as 



< v*v j+k >= f(k/o 



(10) 



Here a\ 
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are the usual Pauli matrices. 



A formal solution for the Green function of Eq. [8] has been derived by Abrikosov and 
RyzhkirJH; their solution generalizes immediately to the commensurate case. The solution 
is expressed in terms of the S— matrix defined for z > z' by 



S(z, z'\ uj) = i^-T z exp I - J dyA(y) 



Here T z is the "space ordering" symbol and 
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(11) 



A = — L> + v ■ a] 

v F 



(12) 



Abrikosov and Ryzhkin show that e.g. the right-moving Green function, Grr, is given for 
z > z' by (we usually do not make the ^-dependence explicit) 



u > : G RR (z, z') = S 22 (oo, z)S 22 (z', -oo)/iS 22 (oo, -oo) 
u < : G RR (z, z') = S 12 (oo, x)S 2 i(z', -oa)/iS u (oa, -oo) 

Further, in this model the density of states for a system of length L is given by 
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(13) 
(14) 



N(u) = -^-TrImlnS{L/2, -L/2; u + it) (15) 



Consider the matrix A (z). Its eigenvalues are 
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±k(z) = ±— J\v{z)\ 2 -u 2 (16) 

Vp 

and so we have 

A(z) = K(z)Q(z)a 3 Q-\z) (17) 

with Q a rotation matrix. 

If k~ x is small compared to the length £ over which A varies, S may be found via a WKB 
approximation. Details are given in the Appendix. The essential result is that S = QRD 
where Q(z) is the rotation matrix which diagonalizes A(z), R is a rotation matrix which is 
close to unity and D is a diagonal matrix which, up to a phase, is 



D = exp -cr 3 J dy 



(18) 



Here 9 = tan~ l {\v\ / uj) and <p = arctan(|i;i|/|v2|). 

From Eq. |15| it is clear that G if k is purely real the density of states vanishes; thus 
obtaining real excitations of energy u requires the existence of regions in which the local 
gap is larger than u. In particular, in models with only phase fluctuations, the local gap 
amplitude \v \ is constant and for |a>| < \v\ the density of states vanishes. This result is of 
course familiar in the context of s-wave superconductivity, where low-lying states and indeed 
dissipation come only from vortices and phase-slip centers, at which the amplitude of the 
superconducting order parameter vanishes. The methods used by Lee, Rice and Anderson!, 
by Sadovskiii, and by many subsequent workers assume that v is Gaussian and peaked at 
\v | = and therefore cannot address the rare amplitude fluctuations which are the relevant 
physics for realistic systems as already noted byi. 

In the limit of infinite correlation length the density of states follows immediately from 
the leading order terms; one finds for a fixed realization of disorder 
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The average density of states is then obtained by averaging this expression over the prob- 
ability distribution of V, leading to the £ — > oo results previously obtained by Sadovskii!. For 
example, for the incommensurate problem the two independently fluctuating components of 
the potential lead to 

N 'Jo Ag 6 v^A^ {M) 

For uj <C Ao, Eq. implies N(u)/N = The analogous expression for the com- 
mensurate problem leads to a density of states proportional to to. 

The expressions presented above suggests that the 'infinite correlation length' results fail 
when ({d z 9) 2 / '«) ~ k, i.e. n^^pp ~ ^ ■ If one considers this estimate as a function of 
frequency, then one sees that in the regions of small potential where low lying states occur, 
k ~ cu. Further, the Lorentzian disorder assumed above has a divergent second moment, so 
{{d z V) 2 ) ~ Aq/(£o) where a is an ultraviolet cutoff of order a lattice constant. Combining 
these factors suggests that the WKB approximation breaks down for u A ~ D 2 (^a). Matching 
this scale to the low energy density of states ~ u 2 yields a residual density of states of order 
£ -1 / 2 . As we shall see in the next section, the numerics suggests rather £ _1 / 3 . Similar 
considerations for the commensurate case would yield N(cu) ~ £ -1 ^ 4 ; as we shall see, the 
numerical results roughly agree. Our arguments suggest that the value of the low energy 
density of states is controlled by fluctuations in the derivative of the random potential. These 
are ultraviolet divergent for the widely used Lorentzian form of the potential fluctuations. 
Study of random potentials with finite second moments and hence derivative fluctuations 
on the scale of £ -1 would be of interest. 

We now consider the density of states at very low energies. Very low energies correspond 
to very long length scales, and at length scales greater than £ one expects the problem to 
map on to one with point-like (delta-correlated) disorder. Ovchinnikov and Erikhmanlii have 
shown that in a fluctuating gap model very similar to the £ — > limit of the commensurate 



potential potential case the density of states diverges as 1/w (with logarithmic corrections 
so the integral is not divergent) and therefore we expect such a divergence in the present 
commensurate potential case also. The interesting question is the dependence on £ of the 
coefficient of the divergent term on the correlation length: this is equivalent to the question 
of the frequency range over which the divergence is visible above the background. 

We present here a qualitative argument indicating that the correct scale is £ _1//2 . We 
note first that the density of states peak in the fluctuating gap model may be traced back 
to the Su-Schrieffer-Heege argument that in such models a change in sign of the back- 
scattering potential produces a mid-gap state. The mean distance between sign changes 
is £ x / 2 . Further, mid-gap states decay exponentially on a scale set by the mean gap so 
the hybridization between mid-gap states may be neglected and they may be treated as 
independent. We conclude that the number of such states per unit length is of order £ -1 / 2 , 
therefore the divergence (which is integrable) must exist for uj < £ -1 ^ 2 . 

This argument may be sharpened. As noted by Bartosch and Kopietz0, at uj = the 
Green function for the commensurate case may be explicitly computed. In the context of 
the WKB formula this may be easily seen from Eq. |TB|: in the commensurate case the fact 
that the potential is purely real means d z ip = while d z 6 — > as uj — > 0, so corrections to 
the WKB result vanish. The number of states at uj = may then be explicitly computed 
from Eq. [15|. At uj = 0, S is purely real except when the potential crosses zero, at which 
point an extra phase m is incurred. Thus in a system of length L the number of states at 
uj = is given precisely by the number of zero-crossings, which is of order L/{a^) 1 / 2 . To 
convert this into an estimate for the density of states as a function of frequency, an estimate 
of the number of nearby states is required. We argue that this may be obtained by scaling, 
assuming uj ~ 1/L; this leads to a divergence ~ l/a;(£a) 1 / 2 . Of course, rough estimates such 
as these will not correctly capture logarithmic terms. 

Bartosch and KopietzS have analyzed the formal solution to the commensurate problem 
in a different way, obtaining an expression which they interpret at the zero-frequency density 
of states per unit length per unit frequency. Their expression diverges exponentially in the 
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size of the system and is consistent neither with the arguments given above nor with the 
numerics to be presented in the next section. In fact their expression closely resembles 
the expression for the real part of Grr which follows from our analysis. We have argued 



states, so we do not consider their results further. 

We now briefly discuss the uo — > density of states for incommensurate potentials. The 
arguments presented above suggest that this will behave differently than in the commensu- 
rate case -the fact that the fluctuating gap has two components means that the zero crossing 
argument is not relevant and the probability for the root mean square gap to vanish is negli- 
gible. We therefore believe that the incommensurate model has a vanishing density of states 
precisely at zero frequency. 



In this section we present results of a detailed numerical study of the density of states of 
the two models. We proceed by writing the real-space form of the Hamiltonian as a matrix, 
choosing a particular realization of the backscattering potential from the distribution defined 
in the previous section, and numerically diagonalizing the Hamiltonian. We obtain the 
density of states by averaging the eigenvalues over an appropriate energy window (of order 
(band width) /50) and then average over many (typically 1000) realizations of the disorder. 
We find that the average over realizations of the disorder converges more rapidly for longer 
correlation lengths than for shorter ones, and converges more slowly for the commensurate 
model than for the incommensurate one. 

We first consider the commensurate model, Eq. [j], which may be diagonalized as it 
stands. For the disorder we chose a Gaussian distribution in which 




that what they have computed is a wavefunction amplitude, not a density of 



III. NUMERICAL RESULTS 



V 3 = Re 5>(g)e ! 



(21) 



and the distribution of V(q) is determined by the kernel 
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K, 



(?) 



A*sinh(l/0 



(22) 



comm 



cosh{l/£) — cos(q — ir) 



These choices correspond to periodic boundary conditions for the fluctuating potential 
and open boundary conditions for the electrons. 

Because H for the commensurate case may be written as a tridiagonal matrix, Sturm- 



(up to L ~ 10 7 ); these system sizes are large enough that boundary effects are entirely 
negligible; study of the length dependence of results has has allowed us to verify that even 
for the smaller systems (L ~ 10 4 ) accessible to direct diagonalization, boundary effects are 
negligible. 

Figure 1 shows the density of states for the commensurate model for several correlation 
lengths and A = 0.2. Several features are immediately evident. First, as noted by other 
workersi, correlation lengths larger than £a = ^f/A are required in order to obtain an 
appreciable pseudogap. Second, the density of states drops only slowly as £ is creased, 
and is surprisingly large even at £ ~ 100. Third, at low energies the density of states is 
approximately constant except for a 'central peak' which is centered at u = and has a 
width which diminishes as the correlation length increases. We have obtained the residual 
(i.e. without central peak) density of states by smoothly extrapolating the calculated N(u) 
to zero, neglecting the upturn. The procedure suffers from ambiguities because our data 
are relatively coarsely binned in frequency and there is some uncertainty about the proper 
functional form for the central peak. The results are shown as filled circles in Fig. 2; the 
dashed line clearly shows that the residual density of states vanishes as with 77 ~ 1/2. 

We now turn to the central peak. It is evident that the width decreases as £ increases. Our 
data are not sufficiently accurate to allow us to determine the precise scaling of the central 
peak with £ and u, but we believe they are consistent with the form N(u, £) ~ f{uj^ 1 ^ 2 ) with 
f(x) given by the Ovchinnikov-Erikhman formEl, f(x) ~ l/xln 3 (x) for x < 1 and / — > 
for x — > 1. In particular, N(u) appears to increase roughly as 1/w and as shown in Fig. 3, 
the area under the central peak scales approximately as £ -1 ^ 2 , although one sees that the 



chain techniquesE-3 may be used to obtain the density of states for very large syst 




em sizes 



10 



corrections to this scaling become appreciable for £ < 100. 

The incommensurate model is defined via a continuum equation; a discretization is there- 
fore necessary. For computational convenience we have adopted 

H = ~ 1 { [ d ]+i,R d j,R ~ ~ M] + V i (4,L d j,R + d ),R d i,i) } ( 23 ) 

3 

with complex Vj gaussian distributed with correlator 

< V;V j+k >= Afc-I*!/* (24) 

Our results for the density of states are shown in Fig. 4. The more rapid drop of N with 
u, expected from the WKB argument, is evident. There is no central peak. Indeed there 
are some indications of a 'central dip' but it is difficult to resolve this question numerically. 
We have obtained the residual density of states by smoothly extrapolating the calculated 
N(u) to zero, neglecting any possible downturn; the procedure is less ambiguous than in 
the commensurate case. The results are shown as filled squares in Fig. Q the scaling with 
£ seems to be closer to £~ 2 / 3 than to the theoretically predicted £ -1 ^ 2 . 

IV. COMPARISON TO APPROXIMATE CALCULATIONS 

In this section we discuss the relation of the numerical results to various approximate 
calculations, in order to obtain insight into the strengths and weakness of the different 
approximations. We begin with the WKB approximation, shown as heavy dashed line in 
Figs. 1 and 4. This is seen to be a good approximation to the calculated density of states for 
not too low energies and not too short correlation lengths; essentially the numerical results 
follows the WKB one until the density of states drops to the residual level shown in Fig. 2. 

We now turn to the continued fraction method of Sadovskii, which for the incommen- 
surate case is compared to numerics in Fig. 4. The qualitative correspondence is seen to 
be good, and to improve for longer correlation lengths. This suggest that the terms which 
Tchernyshyov0 has noted are neglected in Sadovskii's approach are not quantitatively impor- 
tant and become less significant as £ is increased. For infinite correlation length Sadovskii's 
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results are justified by the WKB arguments of section II. Closer examination however shows 
that the low energy large £ behavior of the density of states is not so well represented as one 
can see in Fig. 2, the magnitude differs from the numerical one by factors of order two and 
the scaling with £ is incorrect, being £~ 1 / 2 instead of the numerically determined exponent 
—2/3 to —1. Similar discrepancies arise in the commensurate case, where the Sadovskii 
iV(0) ~ £ -1 / 3 instead of the correct £ -1 ^ 2 . 

We now turn to the two other widely used approximations with a more transparent 
physical content, namely the Lee- Rice- Anderson (LRA) and "fluctuation-exchange" (FLEX). 
The former authors argued that one should approximate the electron self energy by the 
leading order graphs which correspond to the expression 



where Q = 2kp, G = (iu — e p ) 1 is the bare Green's function and K(q) is defined in Eqs. 



for relatively small £ and to a low energy density of states which varies as l/£ for both 
commensurate and incommensurate cases. These incorrect results arise because the self 
energy is too singular; this feature in turn arises because the bare electron Green's function 
is used to describe the intermediate state. 

An alternative scheme is the FLEX method. This is complicated in general but in 
the case of present interest is equivalent to making the Lee-Rice- Anderson calculation self- 
consistent, by using a fully dressed Green's function to compute the one-loop self energy i.e. 
T, FLEX {k,uj) = J G K with G = (iu - e k - T, FLEX {kuj)y l . 

We have numerically solved the FLEX equations for A = 0.2 and various £, the results 
are shown in Fig. |5|. It is clear from the results that the dressed Green function leads 
to a less singular integral and therefore to a larger low energy density of states with less 
£ dependence. These are successes; on the other hand, as with the LRA approach, the 
difference between commensurate and incommensurate cases is lost and the central peak is 
absent. Also as seen in Fig. |5] the location of the peak in N(E) has a strong correlation 




(25) 



|2T| and |22|. This approximation leads to a gap structure which becomes very sharp even 
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length dependence, inconsistent with the numerics and with physical intuition. The value of 
N(0) is found to be larger than that found numerically but the £ dependence is qualitatively 
reasonable (~ £ -1 / 2 )- 

V. CONCLUSION 

To summarize, we have used numerical methods and a WKB analysis of a formal solu- 
tion of a Schroedinger equation to obtain an expression for the Green function of a model 
of a one-dimensional charge density wave in its fluctuation regime. We found that low-lying 
density of states comes from regions where the amplitude of the CDW gap vanishes and we 
emphasize that a proper treatment of a physically relevant model requires a correct treat- 
ment of amplitude fluctuations, which are typically described by a non-Gaussian probability 
distribution, which is difficult to handle either analytically or numerically. 

On the qualitative level, we found that pseudogaps require relatively extreme conditions: 
a drop in the density of states does not begin to appear until the correlation length is 
larger than the basic coherence length i>f/Ao defined by the electron velocity and mean 
field fluctuation amplitude, and the low energy density of states decreases only slowly as £ 
is increased beyond this scale. Physics we have omitted from our model, including quantum 
fluctuations of the pairing potential and the phase space effects characteristic of dimensions 
greater than one, only weakens the tendency to form a gap. The 'pseudogap' observed in 
underdoped high-T c superconductors involves a significant suppression of the low energy 
density of states and therefore implies, at least for these materials, the existence of well 
established, reasonably long ranged pairing fluctuations. 

On the technical side, we have shown that the WKB method (which we suspect can be 
generalized to higher dimensions) and the Sadovskii approximation (which probably cannot) 
provide relatively reasonable estimates of the density of states; other approximations do 
rather poorly, which is unfortunate because they easily generalize well to dimensions larger 
than one. 



13 



Three extensions of this work would be desirable. One is to calculate the conductivity; 
another is to numerically investigate the crossover between the large-£ pseudogap behavior 
and the small-£ constant density of states behavior; the third is to treat non-Gaussian 
amplitude fluctuations. 
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APPENDIX A: WKB APPROXIMATION FOR S 

For z > z', S obeys the equation 

[«9, + A(^)]S = (Al) 

with 

A(z) = iua 3 + A(z) cos (p(z)a\ + A(z) sin p(z)a 2 (A2) 

Here 01,2,3 are the usual Pauli matrices and A and ip are related to the quantities v\ and 
v 2 defined in the text by Vi = A cos ip; v 2 = A sin ip. The eigenvalues of A are ±k(z) with 
k 2 (z) = A 2 (z) — uo 2 If the scale £ over which A varies is much larger than k" 1 , Eq. |A1| may 
be solved via a WKB approximation. Write 

S(z, z') = Q(z)R(z)B(z)I(z') (A3) 

Here Q exp {i6n ■ 8) with 

tan(*0) = A/co (A4) 

ft = (ni, n 2 , n$) = (sin <p, — cosy, 0) (A5) 

and 

D = exp ^{[ K (y) + d 3 (y)]a3 + rfo(y)}j (A6) 

is a diagonal matrix and d 3 and d are functions to be determined. R is a rotation matrix 
which is close to the unit matrix and I expresses the initial conditions. Using Eq. |A3| and 
Eq. RE in Eq. DO gives 



[« + c^RasR" 1 + do - (<9 Z R)RT 1 = [Q _1 0*Q] + k<t 3 (A7) 
Explicitly, [Q _1 <9 2 Q] = i<f • 8 with 
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q = nd z 9 + - sin 29 (ft x z) d z (p + sin 2 9z<p z 



(A8) 



By assuming R is the unit matrix 1 plus small corrections and iterating the equation one 
obtains 



<i 3 = i sin 2 9 d z (p — — 



(d z 9f + -sin 2 (29)(d zV f 



2k I 

R=1 + i(z x q) _ (f x q) 2 + izxd z 



2k 



8k 2 



K 



z x q 
2k 



+ . . 
+ . . . 



(A9) 
(A10) 



This solution may obviously be extended. A non-zero value of do occurs in the third order. 
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APPENDIX: FIGURE CAPTIONS 



Fig. 1 Energy (E) dependence of density of states N calculated numerically for commensurate 
potential model defined by Figs. 1-3 of the text, for different correlation lengths £. The 
mean gap value Ao = 0.2 is shown by a dotted line. The data come from eigenvalues 
binned into intervals of width 0.01. The dashed line shows the 'infinite correlation 
length' result from the WKB approximation. 

Fig. 2 Residual (E — > 0) density of states N calculated numerically for the commensurate 
potential defined by Eqs. 1-3 (filled circles), the incommensurate potential defined 
by Eqs. |10| (diamonds), and the approximation of Sadovskii for commensurate (tri- 
angle) and incommensurate (square) CclSGS, £1X6 plotted against correlation length £ 
on doubly logarithmic axes. Also shown are dotted, dashed, and dash-dotted lines 
indicating scaling with £ -1 ^ 3 , £ -1 ^ 2 , £~ 2 ^ 3 , respectively. For the numerically calculated 
commensurate potential case, the central peak evident in Fig. 1 has been subtracted 
as described in the text. 

Fig. 3 Spectral weight, i.e. integrated area under central peak for the commensurate potential 
case, plotted against square root of inverse correlation length. 

Fig. 4 Energy (E) dependence of density of states N calculated numerically for incommen- 
surate potential model defined by Eqs. 1, 22, 23 of the text, for different correlation 
lengths £. The mean gap value A = 0.2 is shown by the vertical dotted line. The 
results are obtained from numerically calculated eigenvalues binned into intervals of 
width 0.01. The dashed line is the result of the the 'infinite correlation length' WKB 
calculation. 

Fig. 5 Energy (E) dependence of density of states N calculated via the FLEX approximation 
from Eqs. ^ of the text for different correlation lengths and the incommensurate 
potential. The infinite correlation length WKB result is shown by the dashed line. 
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Note that the mean gap value A was fixed at 0.2; the variation of the position of the 
density of states peak is an artifact of the FLEX approximation. 



20 




FIG. 1. 
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